Electric field of a point-like charge in a strong 
magnetic field 



A.E. Shabad* and V.V.Usovt 

*P.N. Lebedev Physics Institute, Moscow, Russia 

E-mail: shabadOlpi . ru 

* Center for Astrophysics, Weizmann Institute of Science, Rehovot, Israel 

E-mail: fnusovOwicc .weizmann. ac . il 

Abstract. We describe the potential produced by a point electric charge placed into 
a constant magnetic field, so strong that the electron Larmour length is much shorter 
than its Compton length. The standard Coulomb law is modified due to the vacuum 
polarization by the external magnetic field. Only mode-2 photons mediate the static 
interaction. The corresponding vacuum polarization component, taken in the one-loop 
approximation, grows linearly with the magnetic field. Thanks to this fact a scaling 
regime occurs in the limit of infinite magnetic field, where the potential is determined 
by a universal function, independent the magnetic field. The scaling regime implies a 
short-range character of interaction in the Larmour scale, expressed as a Yukawa law. 
On the contrary, the electromagnetic interaction regains its long-range character in a 
larger scale, characterized by the Compton length. In this scale the tail of the Yukawa 
potential follows an anisotropic Coulomb law: it decreases as the distance from the 
charge increases, slower along the magnetic field and faster across. The equipotential 
surface is an ellipsoid stretched along the magnetic field. As a whole, the modified 
Coulomb potential is a narrower function than the standard Coulomb function, the 
narrower the stronger the field. The singular behavior in the vicinity of the charge 
remains unsuppressed by the magnetic field. These results may be useful for studying 
atomic spectra in super-strong magnetic fields of several Schwinger's characteristic 
values. 



1. Introduction 

The fact that the vacuum, in which an external magnetic field is present, is an optically 
anisotropic medium, has been known, perhaps, since the time when the nonlinearity of 
quantum electrodynamics was first recognized: in a nonlinear theory electromagnetic 
fields do interact with one another, provided the strength of at least one of them is 
sufficiently strong (of the order of or larger than m 2 /e, where m and e are electron 
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mass and charge, respectively.) If the external field is strong, other fields interact with 
it, the result of the interaction depending upon the direction specified by the external 
field, hence the anisotropy. Depending on the wave amplitude, the electromagnetic 
waves propagation in this medium may be considered as a nonlinear process pQ, 
including the transformation of one photon into two [2] or more photons, or taken 
in the linear approximation with respect to the amplitude. In the latter case, the 
second-rank polarization tensor is responsible for the properties of the medium. In 
the kinematic domain where the photon absorption processes like electron-positron pair 
creation are not allowed, the polarization tensor is symmetric and real, and the medium 
is transparent and birefringent^. In the absorption domain the medium is dichroic [I]. 
The limit of low frequency and momentum belongs to the transparency domain and 
corresponds to a constant anisotropic dielectric permeability of the medium. In this 
limit the polarization tensor may be obtained by differentiations with respect to the 
fields of an effective Lagrangian, calculated on the class of constant external electric 
and magnetic fields. For small values of these fields [H] and for extremely large [I] 
fields the polarization operator was in this way considered using the Heisenberg-Euler 
effective Lagrangian calculated [H] within the one-loop approximation. The knowledge 
of this limit is useful for studying the dielectric screening of the fields that are (almost) 
static and (almost) constant in space. For more general purposes, however, this limit 
is not sufficient, and one should calculate the polarization tensor directly, using the 
Feynman diagram technique of the Furry picture in the external magnetic field. On the 
photon mass shell, i.e. when the photon energy, ho, and 3-momentum, k, are related 
by the free vacuum dispersion law k% = k 2 , such calculations were done by S.Adler 
[H and D.H. Const ant inescu jH]. The results obtained are appropriate for handling the 
photon propagation in weakly dispersive medium, when the dispersion law does not 
essentially deviate from its vacuum shape. The polarization operator for the case of 
general relation between the photon mass and momentum was calculated by I.A.Batalin 
and A.E.Shabad 0, W.Tsai PU, V.N.Baier et al. [II], D.V.Melrose and R.J.Stoneham 
[12]. This gave the possibility of studying the photon propagation [13] under the 
conditions where the deviation from the vacuum dispersion law may be very strong 
either due to the phenomenon of the cyclotron resonance in the vacuum polarization 
[T%] - this phenomenon is responsible for the effect of photon capture by a magnetic 
field [13], [IH], [H] - or due to magnetic fields, much larger than the characteristic value 
B = m 2 /e ~ 4.4 x 10 13 G [IE], [H, [20], [21], or due to the both circumstances [22]. 

Now we start an investigation of problems of electro- and magneto-statics in the 
presence of a strong external magnetic field in the vacuum. For this purpose expressions 
for the dielectric permeability of the magnetized vacuum obtained from the Heisenberg- 
Euler Lagrangian are applicable only as far as the fields slowly varying in the space 

| If the scheme is extended to electron-positron plasma in a magnetic field 3 , the polarization tensor 
is no longer symmetric, but remains Hermitian in the transparency domain. The absorption mechanism 
in this case, apart from the pair creation, includes also the inverse Cherenkov and inverse cyclotron 
radiation of the plasma electrons. 
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are concerned. Otherwise, the spatial dispersion becomes important. For this reason, 
when considering the electric field produced by a point-like electric charge in the present 
paper, we refer again to the polarization tensor calculated off-shell in f|H])-([T%j). taking 
it in the static limit k = 0, but keeping the dependence on k. 

Using the tensor decomposition of the polarization operator and the photon Green 
function, established in 0, [IB] in an approximation- independent way, we find that 
photons of only one polarization mode (mode-2 in the nomenclature of these references) 
may be carriers of electrostatic force. This is in agreement with the fact that the 
electromagnetic field of these photons is, in the static limit, purely electric and 
longitudinal. The photons of the other two modes mediate in this limit the magneto- 
static field of constant currents. 

To describe the static field, produced in the magnetized vacuum by a point electric 
charge, we confine ourselves to considering the asymptotically large magnetic fields 
B ^> B , such that the electron Larmour length (eB)^ 1 ^ 2 is much less than the 
electron Compton length m -1 . Therefore, two different scales occur in the problem: 
the Larmour scale and the Compton scale. A simplifying expression for the mode- 
2 eigenvalue of the polarization operator is used, valid for such fields. It was first 
obtained by Yu.M.Loskutov and V.V.Skobelev [23] within a special two-dimensional 
technique intended for large fields, and by A.E.Shabad [TH], and D.B.Melrose and 
R.J.Stoneham [|l2j as the asymptote of the mode-2 eigenvalue calculated (|H])-([T3]) 
in the one-loop approximation. (See [22] for the detailed derivation of the large-field 
asymptotic behavior.) The most important, now well known, fact about this asymptotic 
behavior (see, e.g., the monographs [21J, [H], [25J) is that the mode-2 eigenvalue contains 
a term linearly growing with the magnetic field. It is sometimes expected that this term 
- it appears in the denominator of the photon propagator and hence of the expression 
for the potential - should lead to suppression of the interaction mediated by mode-2 
photons. In a different problem we already had an opportunity to show that this is not 
exactly the case j2E] , |2Z] • We demonstrate in the present paper that this term is crucial 
for the possibility that a scaling regime might be achieved in the limit of an infinite 
magnetic field, characterized by a comparatively simple universal function, independent 
of the magnetic field. This function gives the potential of a point charge in the energy 
units of (eB) 1 / 2 as a dimensionless function of the space coordinates in the units of the 
electron Larmour length (e-B) -1 / 2 . Except for the closest vicinity of the source charge, 
its form coincides with the Yukawa law (see Eq. ()28)1 below) with the dimensionless mass 
parameter 2a /it (which is the "effective photon mass" discussed in [2E], measured in 
the units (eB) 1 ^ 2 ). Thus, this mass governs the exponential decrease of the potential 
away from its source. As one moves farther from the source, the Yukawa decrease 
ceases, and the potential tail approaches the anisotropic Coulomb shape of the form 
of Eq. ljHKj) . The law of decreasing along the magnetic field is unaffected by the latter, 
the decrease is the fastest in the direction orthogonal to the magnetic field. Therefore, 
the linearly growing term in the mode-2 eigenvalue of the polarization operator leads, 
first, to the faster decrease of the potential in the direction across the magnetic field 
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for large distances, and, second, to its steeper shape for small distances. This may be 
recognized as suppression, indeed. On the other hand, the long-distance behavior along 
the magnetic field, as well as the standard Coulomb singularity near the source do not 
sense the magnetic field at all, no matter how strong it is. 

In Section 2 we give an approximation-independent derivation for the form of 
the potential of a point charge in terms of the mode-2 component of the photon 
propagator. In Section 3 we consider the short-distance behavior of the potential. 
Three terms of its asymptotic expansion near the source are explicitly written. The 
universal function for the scaling infinite-field limit is obtained as a one-fold integral. 
In Section 4 the anisotropic Coulomb law is obtained for large distances by considering 
different mathematical means, appropriate for different remote regions in the space. The 
analytical results for the potential are supplemented by its shapes drawn by a computer, 
presented in Figures 1-5. 



2. General representation for the static potential of a point-like charge 

Electromagnetic 4- vector potential produced by the 4-current j v {y) is 

A,(x) = J D^(x - y)f(y)d 4 y, fi,,u = 0, 1, 2, 3. (1) 

Here x and y are 4-coordinates, and D^ix — y) is the photon Green function in a 
magnetic field in the coordinate representation. The metric in the Minkowski space is 
defined so that diag = (1, —1, —1, —1). Eq. ((TJ) defines the 4-vector potential with 
the arbitrariness of a free-field solution, including the gauge arbitrariness. If the photon 
Green function is chosen as causal, only the gauge arbitrariness remains. 
The current of a point-like static charge q, placed in the point y = is 

/(!/)=?W 3 (y). (2) 

Hence 

/oo 
£Vo(zo - y ,x.)dy 
-oo 

/OO POO 
D^{x Q + y , x)d?/ = q D^(y , x.)dy . (3) 
-OO J —oo 

This 4-vector potential is also static. 

Where there is no magnetic field, and the photon propagator is free and taken in 
the Feynman gauge (with the pole handled in a causal way) 

D^x -y) = D%{x -y) = i4 ^_ y)2 , (4) 
only the zeroth component of the 4-vector potential is present: 

»?(*) = , £ D^ B , x)* ^£^4^,. (5) 
This is the Coulomb potential in the Heaviside-Lorentz system of units used throughout. 



Electric field of a point-like charge in a strong magnetic field 



5 



Let there be an external magnetic field B directed along axis 3 in the Lorentz frame 
where the charge q is at rest in the origin x = , and no external electric field exists in 
this frame. Call this frame special. Define the Fourier transform as 

' _ I nvnl, / ■ r \ I , x 



exp(ikx)D IMU (k) d 4 k, 



0,1,2,3. 



(6) 



Then (J2J) becomes 

q 



(27T) 



exp(i(k y -k^))D^(k) d 4 My = j D^O, k) exp(-ikx)d 3 A;. (7) 



The four 4-eigenvectors b^\ a = 1,2, 3, 4, of the polarization tensor H"^ (HI HBl IHj 
take in the special frame (and arbitrary normalization) the form - the components are 
counted downwards as u = 0, 1, 2, 3 - 

h 

/ k \ 



6« 



/ 



\ 



( \ 

V / 
o \ 

k 2 

-h 





+ ki 






\k J 



(8) 



h 

k 2 

V ** / 

Among them there is only one, whose zeroth component survives the substitution ko 

( 2 ) 

It is b v . This implies that out of the four ingredients of the general decomposition of 
the photon propagator 



0. 



&£° b^ 



D a {k) 



a=l 

-(e + K a (k)y\ 

arbitrary, 



a- 



=1,2,3 



(9) 



where n a (k) are scalar eigenvalues of the polarization tensor: 

n;(k)b^ = Ka (k)bp, K 4 (k) = o, (io) 

only the term with a = 2, D 2 (k)b^b^ /(b^) 2 , participates in ((Zj), i.e. only mode- 
2 (virtual) photon may be a carrier of electro-static interaction, and not photons of 
modes 1,2, nor the purely gauge mode 4. Bearing in mind that (b^) 2 = k\ — k 2 ., we 
have 



A)(x) 



Q 



(27T) 



-ikx j3 



d 6 k 



^i,2, 3 (x) = 0. 



k 2 -K 2 (0,k 2 ,kl) 

Here k\ = k\ + k\. Thus, the static charge gives rise to electric field only, as it 
might be expected. The gauge arbitrariness in the choice of the photon propagator 
D±(h)bpbv = D^tykpkv indicated in (JUJ) has no effect in (JTTj). Certainly, the potential 
(fTT|) is defined up to gauge transformations. 
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The result that only mode-2 photons mediate electrostatic interaction may be 
understood, if we inspect electric and magnetic components of the fields of the eigen- 
modes obtained from their 4-vector potentials (jHJ) in the standard way: e^ = koh^ a ' — 
kb a) , b>> = k x b^. These are [T3| 

e (1) = -^Ab, hW = x k 3 ), (12) 

k± k± 

ef = k ± k 3 , ef = |(fc 3 2 - kl), = - ko (k ± x |) , (13) 

e< 8 > = -ko ( ^ x , h± W = -£:*,, h3® = ^ ^. (14) 

Here the cross stands for the vector product, and the boldfaced letters with subscripts 
3 and _L denote vectors along the directions, parallel and perpendicular to the external 
magnetic field, resp. 

The photon energy and momenta here are not, generally, connected by any 
dispersion law. Therefore, we can discuss polarizations of virtual photons - carriers 
of the interaction - basing on Eqs. (|T2"|) - (jl4|) . The electric field e in mode 1 is parallel to 
kx, in mode 2 it lies in the plane containing the vectors k, B, in mode 3 it is orthogonal 
to this plane, i.e. mode 3 is always transversely polarized, e^k = 0. For the special case 
of the virtual photon propagation transverse to the external magnetic field, k 3 = 0, (this 
reduces to the general case of propagation under any angle 6 ^ by a Lorentz boost 
along the external magnetic field), mode 2 is transversely polarized , e^k = 0, as is 
always the case for mode 3. Mode 1 for transverse propagation, k 3 = 0, is longitudinally 
polarized, e^ 3 -* x k = 0, and its magnetic field is zero. The lowest-lying cyclotron 
resonance of the vacuum polarization ^3], the one that corresponds to the threshold 
A)q — fc| = 4m 2 of creation of the pair of electron and positron in the lowest Landau 
state each, belongs to mode 2. It gives rise to the photon capture effect with the photon 
turning into a free jEj or bound fHl El electron-positron pair. Another consequence 
of the cyclotron resonance is that a real photon of mode 2 undergoes the strongest 
refraction in the large magnetic field limit [22] even if its frequency is far beyond the 
pair production threshold. 

In the static limit ko = the magnetic field in mode 2 disappears, h.^ = 0, while 
its electric field is collinear with k, e^ 2 ^ = k. It becomes a purely longitudinal virtual 
photon. Unlike mode 2, in modes 1 and 3 in the static limit ko = only the magnetic 
fields survive: e^ 1 ' 2 ) = 0, h« = k ± x B, hf = -k ± k 3 , h { 3 3) = k 2 ± , h^k = 0. (Here 
normalizations are arbitrary and kept fixed only within the same mode). A virtual 
mode-1 photon carries magneto-static interaction. It is responsible for magnetic field 
produced by a current flowing through a straight-linear conductor oriented along the 
external magnetic field. In accordance with the above formula for h^ 1 ) its magnetic field 
is orthogonal both to B and to the radial direction in the transverse plane k^, along 
which the magnetic field of the current decreases. The mode-3 photon contributes as an 
interaction carrier in the problem of a magneto-static field produced by a solenoid with 
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its axis along axis 3. In the present paper, however, we do not consider magneto-static 
problems. 

In the asymptotic limit of high magnetic field B 3> fcf, B ^> m 2 /e = Bq the 
eigenvalue K 2 (0,k), as calculated within the one- loop approximation [§], JH|, with the 
accuracy of terms that grow with B only as its logarithm and slower is ^H] 

where b = B/B Q and 

T(y) = y f = 1 - -^J_ m (16) 

w y J 1 + 2/(1 - 77 2 ) 2VRT+^) v/T+i-v^ 

Note that T(y -> 0) ~ 2y/3, T(oo) = 1. 

Other eigenvalues, Ki^, do not contain the coefficient b that provides the linear 
increase of f)15jl with the magnetic field. 

Expression for k 2 ([To]) is to be used in (|TT|) or, equivalently, in the expression 

^ = 2^/ ^(L gTlg^TO '^ (17) 

obtained from (|11|) by integration over the angle between the 2- vectors kj_ and Xj_, 
which are projections of k and x onto the plane transverse to the magnetic field. In 
(fTTj) x± = \Jx\ + x\ and J is the Bessel function of order zero. 

The nontrivial - other than the leading asymptote ~ fc| near the point k 3 = k± = 
- dependence of the polarization operator eigenvalue (|TH|) on the photon momentum 
components k$, k± is the spatial dispersion. We shall see in Section 3 that it is important 
in the vicinity of the charge, where the field has a large gradient. As for the anisotropic 
behavior far from the charge, to be studied in Section 4, only the above asymptote is 
essential, inferable also from the Heisenberg-Euler Lagrangian. 

3. Short-distance behavior 

3.1. Asymptotic expansion 

To consider the behavior of the potential near its point-like source let us add to and 
subtract from (fTTjl the standard Coulomb potential (jHJ) in the form 

A c q f e" ikx d 3 A; 



K J ~ (2vr)3 J k 2 



J j (* ± * ± )(£L^)d*» 



4 



so that 



2(2*)* J q ^ ^ Vi-oo kl + kj J ± 4vr 
A (x) = A c (x)-AA (x), (19) 
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where 



AAo(x) 



J (k±x±) / - dfc 3 dA; 2 . (20) 



2(2tt) 2 7 uv ^ x ' 7_ 00 V k\ + fc 3 2 fci + k\ - k 2 (0, A; 2 
Note that the function AAq(x3,x±) is an entire function of x±, since the exponential 
in (|15|) provides convergence of the integral (ffllj) for any complex value of this variable. 
Keeping quadratic terms in the power series expansion of Jo(k±x±) and exp(— ik^x^) in 
(|2()jl we obtain the first three terms of the asymptotic expansion of the potential (j!7j) 
near the origin X3 = x± = 

Ao(x) ~ -p- (A - 2m (^ - (2mx ± ) 2 C ± - (2mx 3 ) 2 C|i) J , (21) 
47r V |x| y 

where C, C±_ and Cy are dimensionless positive constants depending on the external 

field: 

C = — AA (0) = 
qm 



7T 2 







«6m r° T /j_\ r ex pj-^J d ^ 



0± 16mWo U^Vio (*i + + *§ " ^(0, *§, 1 ' 



_ ab r (J4 \ r feaexp^-^jdfc, 

Thanks to the exponential factor the integrals over k\ here are fast converging. The 
resulting functions decrease for large k% as ~ I//C3, so the remaining integrals over k% in 
(|2*2*|) . (J2*3J| . (|2*4*j) converge, bearing in mind that T is a bounded function. The values of 
the coefficients (|22jl. ([23)1. calculated for four values of the magnetic field 6 = 10 4 , 
b = 10 5 , 6 = 10 6 and b = 10 10 are, respectively, C = 2.21, 9.08, 31.37, 32.70 x 10 2 , 
C ± = 75.9, 2.58 x 10 3 , 8.38 x 10 4 , 8.49 x 10 10 , C\\ = 174.3, 5.55 x 10 3 , 1.76 x 10 5 , 
1.67 x 10 11 . 

In Figs. [U and 121 functions A (xs,0) and A (0,x±) are drawn by a computer using 
(fT7j) for three values of the magnetic field: b = 10 4 , b = 10 5 and b = 10 6 . In agreement 
with the above analysis the curves in Figs. ni2 a PP roac h the Coulomb law q/47r\x3\ or 
q/4nx± as X3 — > or x± — > 0, respectively. In Fig. El the ratio of the modified Coulomb 
law (fTTj) to the standard Coulomb law © at x 3 = is plotted against the transversal 
coordinate x± for the magnetic field B = 10 4 i?o- It deviates from unity as we move 
away from the origin x±_ = 0. The dotted line is the interpolation corresponding to 
the potential taken in the Yukawa form (q/Airx ±) exp(—2mCx±) with C defined by Eq. 
(122(1 to be in this case C =2.21. This interpolation is in agreement with the first two 
terms of the asymptotic expansion (|21|). The Yukawa interpolation is good within the 
range < x±_ < 0.1(2m)~ 1 = 5L B . 
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Figure 1. Modified Coulomb potential ^0(2:3, 0) i|17[l plotted along the axis X3 passing 
through the charge q parallel to the magnetic field. Dashed lines correspond to three 
values of the magnetic field (from left to right): B = W 6 B , B = 10 5 B and B = 10 4 B . 
Thick solid line is the standard Coulomb law © Aq(x3,0) = q/Anxs. Dashed lines 
approach asymptotically the solid line at the both edges of the figure. The abscissa 
represents the distance along the magnetic field in the units (2m)" 1 . The ordinate 
represents the potential in the units qm /2tt. 



3.2. Scaling regime 

It is remarkable to note that in the infinite-magnetic-field limit, potential (fT^j) . measured 
in the inverse Larmour length Lg 1 = \J eB units becomes a universal, magnetic-field- 
independent function of coordinates measured in Larmour units Lb- To establish this 
scaling regime let us make the change of variables in the integral ()20|) k 3 = k 3 L-Q, 
k± = k±LB and define the new dimensionless coordinates 23 = x 3 Lb, x± = x±Lb- After 
that, we may substitute unity in place of the function T(fcf/4m 2 L B ), since it is the 
value of this bounded function in the limit (eB/Am 2 ) = l/4m 2 L B = 00. Then Eq. fl2()j) 
becomes 

AA (x) ~ I \ J (k ± x ± ) / fr ^ r - dhdkl. (25) 

2(2tt)2L b J q yk 2 ± + kl p ± + k 2 + ^ exp J 

The /^-integration here can be performed by calculating the residues in the poles on the 
imaginary axis. Finally one gets 
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Figure 2. Modified Coulomb potential Aq(0,x±) (|17|l plotted along the axis x± 
passing through the charge q transversely to the magnetic field. Dashed lines 
correspond to three values of the magnetic field (from left to right): B = 10 6 Bo. B = 
10 5 B and B = 10 4 i?o- Solid line is the standard Coulomb law © A% (0, x±) = q/^-KX^. 
Dashed lines approach asymptotically the solid line at the left edge of the figure and the 
short solid lines A o (0,x±) = q/^nx'^ at the right edge. The abscissa represents 
the distance across the magnetic field in the units (2m) -1 . The ordinate represents 
the potential in the units qm/2n. 




This simple representation can be further simplified if 23 or x± are large in the Larmour 
scale: |x 3 | ^> 1, or \x±\ ^> 1. In this case the integration in (|26|) is restricted to the 
domain fc^ < 1 where the exponential exp(— k\/2) should be taken as unity. Then 
is reduced to 



, . q 1 -exp{- {^Y Jx 2 + x 2 } , 

aa (x) ~ -± P1 — -> ( 2? ) 



which implies the Yukawa law for (fTTJj) 



4,(x) ~ ytf(x) = * lv v L 3J = ^- iV V/ 1 IJ . (28) 

1 ^ 01 ^ 4ttL b y/xl + 3% 4tt |x| 1 ; 

The latter fact can be established if we return to the previous representation ([25)1 . 
which can then be traced back to (|19|) . (|2(Jj) and finally to (|17|) with 

2a 2ab 
-k 2 (0,oo,0) = — T = m 2 = M 2 (29) 

substituted for — /^(O, fcf, kj_) in the denominator. Here M is the "effective photon 
mass" of Ref.[2Hl- The deviation of (|26J) from ()27|) when x± and £3 are both small is 
not important in the sum (j!9)) against the background of the first term Aq (x), singular 
in the origin x± = 23 = 0. Therefore, the Yukawa law ()28|) is fulfilled in the vicinity of 
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Figure 3. Ratio of the modified Coulomb potential l|17fl to the standard Coulomb 
potential JSJ Aq(0, x±)/Aq (0, x±) drawn as a solid line in function of the transverse 
coordinate x± for the magnetic field value B — 10 4 -Bo- The dashed line is the constant 
(x±/x'_ L ) = (f + aB/Z-KBn)- 1 / 2 = 0.338 from the law (jSSJ). The dotted line is the 
Yukawa interpolation. The abscissa is the same as in Fig|3 



the charge. By comparing (|27|) with the second term of the expansion (|2*T|) we deduce 
that for strong fields 



asymptotically. For the four values of of the external field b = 10 4 , 10 5 , 10 6 , 10 10 the 
values of C calculated following flUJ) are: 3.41, 10.78, 34.01 and 34.08 x 10 2 . We face 
the coincidence with its values listed in previous Subsection, the better the higher the 



The Yukawa law (|2*5j) establishes the short-range character of the static 
electromagnetic forces in the Larmour scale. To avoid a misunderstanding, stress that 
the genuine photon mass understood as its rest energy is always exactly equal to zero as a 
consequence of the gauge invariance reflected in the approximation-independent relation 
K a (0, 0, 0) = 0. Correspondingly, the potential, produced by a static charge, should be 
long-range for sufficiently large distances. This is the case, indeed. The point is that the 
scaling regime ([26|1 . as it follows from its derivation, is valid in the limit eB /4m 2 = oo. In 
this limit the ratio of the Compton length (2m) _1 to the Larmour length Lb = (eB)^ 1 ^ 2 
is infinite. Therefore, Eqs. (I2fij) . ()28)1 cannot be extended to the distances of Compton 
scale and larger, where the long-range character of the electromagnetic interaction is 
restored, as we shall see in the next Section. 




(30) 



field. 
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4. Large-distance behavior 

4-1- Large x± in Larmour scale 

For large transverse distances the term linearly growing with the magnetic field (JT5j) 
leads to suppression of the static potential in the transverse direction. 
To be more precise, consider the region 

X± y> Wb = ^2 (31) 
Once the Bessel function J in (JT7|) oscillates and decreases for large values of its 
argument k±x±, the main contribution into the integral over k\ in (jl7j) comes from 
the integration variable domain k\ <C 2m 2 b, and the dependence upon k\ in K2 may 
thus be disregarded. Then the /c^-integration in (fTTjl can be explicitly performed to give 
(we use Eq.6.532.4 of the reference book |29j) 

A {x 3 , x ± ) ~ /C (x ±y Jk% - K 2 (0,fc§,0)^ cos(k 3 x 3 )dk 3 , (32) 

where /Co is the McDonald function of order zero, and 

K2 (0,40) = -^m 2 T(lj. (33) 

As the McDonald function /Co decreases exponentially when its argument increases, 
only small values of the square root contribute into integral (J32|) . and the /^-integration 
domain in it is restricted to the interval k 2 <C 4m 2 , wherein 

(^4m 2 ) 6m 2 ^ ^ 

Then the potential form (|32j) becomes (we use Eq. 6.671.14 of the reference book 29 ) 

2q f°° ( I ab\ 

A (x 3 , x±) ~ -^-^ J /C I x ± k 3] / 1 + — J cos(A;3X3)dA;3 = 

1 Q if, ab\ l/2 . , or , 
= =, x,—x±ll-\ , x, > x±. 35) 

Eq. (j35J) is an anisotropic Coulomb law, according to which the attraction force 
decreases with distance from the source along the transverse direction faster than 
along the magnetic field (remind that b = (B/B ) ^> 1), but remains long-range. 
The equipotential surface is an ellipsoid stretched along the magnetic field. The 
electric field of the charge E = — VA (x 3 , x±) is a vector with the components 
(q/2ir)(xl + /5 2 x 2 L )~ 3 / 2 (x 3 , /? 2 x_l), where (5 = (1 + abj'i-n) 1 ! 2 . It is not directed towards 
the charge, but makes an angle with the radius- vector r, cos0 = {x\ + f3 2 x 2 L )(x 2 + 
P 4 x 2 L )~ 1 ^ 2 (x 2 + x^) -1 / 2 . In the limit of infinite magnetic field, (3 = 00, the electric field 
of the point charge is directed normally to the axis x 3 . This regime corresponds to the 
dielectric permeability of the vacuum independent of the frequency, with its dependence 
on k (spatial dispersion) being reduced solely to that upon the angle in the space (cf, 

[22]). 
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The result (J35|) is in agreement with the curves in Figs|2l 121 drawn for x 3 = 0. The 
curves in Fig|2 approach the standard Coulomb law when x± — > in accordance with 
121)) . as explained above, rather sharply fall down, following the Yukawa law (|2*H|) in the 
Larmour range x± ~ Lb to reach the asymptotic long-range regime A (0,x±) ~ q/Atrx^ 
for larger x± in the region (|31j) . 

4-2. Large x 3 

It remains to consider the remote coordinate region of large x 3 , complementary to ([31)1 . 

To this end we apply the residue method to the inner integral over k 3 in 1)17)) . Using 
the integral representation ()16)) the function n 2 (fTB^) may be, for a fixed positive value of 
kj_, analytically continued from the real values of the variable k 3 into the whole complex 
plane of it, cut along two fragments of the imaginary axis. In the lower half-plane the 
cut runs from Im k 3 = —2m down to Im k 3 = — 00, while in the upper half-plane it 
extends within the limits 2m < Im k 3 < 00. Other singularities of the /^-integrand 
in ()17)) are poles yielded by zeros of the denominator, i.e. solutions of the equation 
(associated with the photon dispersion equation) 



As k± varies within the limits (0, 00) two roots of this equation kf = ±iK(k±) move 
along the imaginary axis from the point K(0) = to the points kf = HK(oo) = ±i2m 
[T3] [TB"] . [22] . There is yet another branch of the solution to equation (pjj). corresponding 
to the photon absorption via the 7 — > e + e _ -decay, but the corresponding poles lie in 
the nonphysical sheet of the described complex plane, behind the cuts, and will not be 
of importance for the consideration below. 

Let us consider positive values of x 3 . Negative values can be handled in an analogous 
way. Turning the positive part of integration path < k 3 < 00 clockwise to the lower 
half-plane by the angle tt/2, and the negative part —00 < k% < counterclockwise by 
the same angle, and referring to the fact that the exponential exp (— i/^a^) decreases, for 
£3 > 0, in the lower half-plane of k% as | A^3 1 — > 00 so that the integrals over the remote 
arcs may be omitted, we get a representation for the inner integral in ()17j) 



where Res(/r[) designates the residue of the expression D2(0, — \ks\ 2 , Ar[) = 
(Jfei- l^p-z^OH^iUi))" 1 in the pole k^ = -iK(k 2 ± ), while A(|A; 3 | 2 , k\ ) = 
L?2 (0, — 1^3 1 2 + i0, kjj — D 2 (0, — \k 3 \ 2 — iO, kjj is the cut discontinuity. It was explained 
above that < K{k 2 1 ) < 2m everywhere but in the limit k±_ —>■ 00, where K = 2m. 
Consequently the residue term in ()37)) dominates over the cut-discontinuity term 
everywhere in the fc^-integration domain in the outer integral in (j!7j) . except for the 
region near k± = 00. In this limit, however, k 2 disappears due to the exponential in 



k 2 ± + k 2 3 -K 2 (0,k 2 3 ,k 2 ± )=0. 



(36) 




(37) 
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(JT5J), together with the cut discontinuity, since the latter is only due to the branching 
points in the function (jl6j) . Therefore, keeping the residue term in ()37|) as the leading 
one, we neglect the contribution that decreases with large longitudinal distance at least 
as fast as exp(— 2m|x 3 |). In this way we come to the asymptotic representation of the 
potential (fTTj) in the region of large longitudinal distances \x\s ^> (2m) _1 (the negative 
23 at this step are also included - to treat them the fragments of integration path were 
rotated in the directions opposite to the above) 

AJ-x) ~ JL r M k ±X±) p ~(K(ki)x 3 ) dk 2 f oox 

A °^-8nJ K(kl)(l + H(-K*(kl),kl) e ^' (38) 

where 

/j9 ,,. 2abm 2 ( k 2 , \ d m ( k\ \ 

Here K 2 (k\) is the solution of equation ()36|) in the negative region of the variable /c| - 
see j22 for its form. K(oo) = 2m, K(0) = 0. T is given by (fToT) . 

Due to the exponential factor in the integrand of (J38j) . for large x% the main 
contribution comes from the integration region of k± that provides minimum to the 
function K(k±). The minimum value of K(k±) is zero. It is achieved in the point 
k± = - a manifestation of the fact that the photon mass defined as its rest energy 
is strictly equal to zero owing to the gauge invariance: K a (ko = k% = k = 0) = 0. 
In view of ffSlIj) and (jM)), near the point k± = the dispersion equation (J35|) has 



the solution K(k±) = k±/^ 1 + ab/3ir. Simultaneously, near the minimum point 
1 + H(0, 0) = 1 + ab/3Tr. With these substitutions and the use of 6.661.1 of [29j, 
Eq. becomes again the anisotropic Coulomb law (q/4ir) /((x' ± ) 2 + X3) 1 / 2 . We 
have, therefore, established its validity everywhere in the region remote from the center, 
irrespectively of the direction. In agreement with this result the curves in Figure 1 for 
^4o(^3,0) are approaching the Coulomb law q/A^x^ as x% grows. To see how fast they 
are doing this, set x± = in the difference ([20)1 . By replacing the exponential in the 
expression for k 2 (fTB^) by unity we overestimate |At4 (x3, 0)|. Then the integration over 
k\ can be explicitly performed to yield the inequality 



|A4,(s 3 ,0)| < 



(40) 



The only singularities of the integrand of ()40)1 as a function of the complex variable k% 
are (inverse square root) branchings on the imaginary axis in the points ±i2m owing 
to the analytical properties of the function T(y) f|16j) . Note that the argument of the 
logarithm in (|4()jl may disappear only on the second sheet of the complex plane. By 
turning the integration path to the lower half-plane, when £3 is positive, or to the upper 
half-plane otherwise, we reduce (|4()|) to an integral of the cut discontinuity multiplied 
by exp(x3lmfc 3 ) along the imaginary axis of k$, the minimum integration value of |Im/c 3 | 
being 2m. Therefore the difference ()20j) between the potential A (x^, 0) and its large-23 
asymptote g/47r|x 3 | decreases in Figure 1 at least as fast as exp(— 2m|x 3 |). 

Eq.(|3*H|) was used for computer calculation with the large value x 3 = 10m _1 . It has 
lead to the curve shown in FigslU El In the region (|3*T]) it agrees with the result (f33J) . 




Figure 4. Modified Coulomb potential l|38fl plotted against the transverse coordinate 
x± with the longitudinal coordinate fixed at the large value x$ = 10(2m) _1 . The 
dashed line corresponds to the magnetic field value B = 10 4 _Bo- Solid line is the 
standard Coulomb law © Aq(x 3 ,x±) = qm/2Tr((2mx±) 2 + lOO) -1 / 2 . The dashed line 
is indistinguishable from the anisotropic Coulomb law Ij35(l in the scale of the drawing. 
The coordinate axes are the same as in Fig|21 
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Figure 5. The same as FigQ] but viewed at a detailed scale near x± — 0. The thin 
solid line is the anisotropic Coulomb law l|35|) (?m/27r((0.338 • 2mxj_) 2 + lOO) -1 / 2 



valid in that region (Lb = 0.02(2m) 1 for b = 10 4 ). In practice and (|H5|) are the 
same. 
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5. Conclusion 

The modification of the Coulomb law should affect, first of all, the field of an atomic 
nucleus, placed in a magnetic field. Properties of matter (including individual atoms) 
at the surface of strongly magnetized neutron stars and various physical processes (such 
as radiation of particles) where the electric field of particles is important (for a review 
on physics of strongly magnetized neutron stars, see [30J ) may become sensitive to the 
present modification of the Coulomb law. 
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